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ABSTRACT 


Dynamic  buckling  loads  are  obtained  for  axlsynmetrlc  spherical  caps  with 
initial  imperfections.  Two  types  of  loading  are  considered,  namely,  step 
loading  with  infinite  duration  and  right  triangular  pulse.  Solutions  of  per- 
fect spherical  caps  under  step  loading  are  in  excellent  agreement  with  previous 
findings.  Results  show  that  initial  imperfections  do  indeed  have  the  effect 
of  reducing  the  buckling  capacity  for  both  dynamic  and  static  responses, 
although  they  are  affected  in  a different  manner.  From  the  solutions  obtained 
for  triangular  pulse  situations,  it  is  reveed.ed  that  pulse  duration  has  a 
very  significant  impact  on  the  magnitude  of  the  dynamic  buckling  load.  When 
comparing  these  solutions  with  those  of  step  loading,  it  is  concluded  that 
the  step  loading  with  infinite  duration  is  the  limiting  case  of  a triangular 
pulse,  and  that  the  step  loading  provides  the  most  severe  loading  situation 


for  dynamic  analysis. 


INTRODUCTION 


Dynamic  buckling  analysis  of  shell  structures  has  received  considerable 
attention  in  the  literature.  Shell  structures  designed  according  to  quas^- 
static  analysis  may  fail  under  conditions  of  dynamic  loading.  It  is  also 
found  that  initial  imperfections  in  spiv  rical  shell  structures  have  the  effect 
of  reducing  shell  static  buckling  capacity  (Refs.  ^l,2,3,i*j).  Initial  imper- 
fection of  shell  structures,  in  fact,  is  sometimes  used  to  account  for  the 
discrepancy  between  results  obtained  from  experimental  tests  and  theorectical 
analyses . 

The  problem  of  axisymmetric  dynamic  snap-through  of  clamped  spherical 
caps  under  implusive  loading  was  first  solved  by  fTUi-.pj.rey  and  Bodner  The 

same  problem  under  instantaneously  applied  step  loading  was  studied  by 
Budiansky  and  Roth  ^6^  and  Simitses  Rayleigh-Ritz  or  Galerkin  methods 

were  used  tor  shell  dynamic  buckling  analysis  in  those  studies.  Archer  and 
Lange  discussed  the  same  problem  numerically  by  solving  governing  differen- 
tial equations  by  a combination  of  finite-difference  technique  and  Potter 
method.  Lock,  Okubo  and  Whittier  approached  this  problem  by  performing 
experiments  for  a particular  shell  geometry. 

The  same  problem  under  a uniform  step  pressure  is  solved  numerically 
in  Refs.  ^io-iuj  and  results  obtained  from  these  investigations  are  in 
reasonable  agreement.  From  this  agreement,  it  may  be  said  that  the  axisym- 
metric dynamic  buckling  criterion  - dynamic  snap-through  - suggested  by 
Budiansky  and  Roth  has  generally  been  accepted. 

Asymmetric  dynamic  buckling  analysis  of  spherically  caps  is  examined  in 
Refs.  [13.15]  in  which  a perturbation  approach  was  utilized  to  deal  with 
asymmetric  deformation  mode. 


* Numbers  in  brackets  designate  references  at  end  of  paper. 


Effect  of  Initial  imperfections  on  static  shell  buckling  analysis  was 
studied  by  Hutchinson  ^1^  who  employed  the  Koiter  approach  to  determine  the 
approximate  asymmetric  buckling  load  for  spheric?"'  shells  with  initial 


imperfections.  Koga  and  Hoff 


H 


gave  the  values  of  buckling  pressure  for 


complete  shells  with  axisymmetric  dimple  type  imperfections.  Numerical  solutions 


for  asymmetric  buckling  loads  are  given  in  Refs. 


for  spherical  caps  with 


initial  asymmetric  Imperfections. 

In  this  paper,  attention  is  focused  on  the  effect  of  initial 
Imperfections  on  the  dynamic  buckling  of  axisymmetric  spherical  caps.  The 
influence  of  different  types  of  dynamic  pulse  shape  on  shell  buckling  capacity 
is  also  considered. 

In  the  next  section,  the  governing  equations  and  solution  methods  utilized 
in  the  present  paper  are  discussed.  This  is  followed  by  a description  of 
dynamic  buckling  criteria.  To  verify  the  present  approach,  a comparison  of 
the  present  results  with  those  in  the  literature  is  given  for  axisymmetric 
buckling  loads  of  spherical  caps  under  uniform  step  pressure  and  with  initial 
imperfections  of  different  magnitude.  Also  presented  are  the  solutions 
associated  with  a triangular  pulse  with  different  time  durations.  Conclusions 
and  a general  discussion  are  given  in  the  final  section. 


GOVERNING  EQUATIONS 


The  gecanetry  of  a clamped  spherical  cap  is  shown  in  Fig.  1,  in  which 


H is  the  central  height  and  R denotes  the  shell  radius;  a is  the  base  radius 


W(r,t)  6ind  U(r,t)  are  displacement  components  along  normal  and  tangential 


is  the  initial  imperfection 


directions  at  time  t 


The  undeformed  shape  of  the  perfect  shell  can  be  adequately  described  by 


where  r is  the  radial  coordinate 


The  general  differential  equations  for  the  response  of  spherical  caps 


with  initial  imperfections  are  given  in  Ref 


only  the  situation  of  the  axisymmetric  deformation  with  the  inclusion  of 


dynamic  effects.  The  equilibrium  equations  associated  with  this  situation  are 


where  7^ ()=()"+()' /r;  D = Eh3/l2(l  - ) , E is  Young's  modulus 


h is  the  shell  thickness,  v is  the  Poisson's  ratio,  and 


unit  volume  of  shell;  prime  and  dot  denote  differentiations  with  respect 


to  r and  t,  respectively,  and  0 is  the  circumferential  coordinate 


the  applied  loading  and  + W;  stress  resultants  are  related  to  strains  by 


in  which  strains  are  expressed  in  terms  of  displacements  by 


The  terns  on  the  right  hand  sides  of  Eqs.  (2)  are  the  inertia  forces 


due  to  normal  and  radial  displacements  of  the  shell.  For  convenience,  the 


following  nondimen sional  quantities  and  operations  are  introduced 


q^j.  = UEh/R^m^ 
p(r,t)  = q(r,t)/q 


where  q^j.  is  the  classical  buckling  pressure  of  a complete  spherical  shell  of 


With  the  adoption  of  Eqs.  (5),  the  nondimensional  forms  of  Eqb . (2) 


Where  e eind  Eq  are  nondimensional  quantities  of  e_  and  Eg,  respectively 


It  is  noted  that,  in  view  of  the  assumed  shallowness,  the  effect  of 
radial  inertia  force  is  neglected  in  Eq.  (6b).  Eqs.  (6)  constitute  the  basic 
equations  for  the  analysis  subjected  to  prescribed  boundary  conditions. 
Boundary  condition  considered  in  this  paper  are  clamped: 

u(1,t)  = w(1,t)  = w'(1,t)  = 0.  Due  -to  symmetry  at  the  apex,  we  also 
have  u(0,t)  = w'(0,t)  = 0. 
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METHOD  OF  SOLUTION 

Eqs.  (6)  are  not  only  lengthy  but  also  highly  nonlinear  in  w,  the  major 
variable.  The  solution  plan  is  to  superpose  a finite-difference  mesh  on  the 
one-dimensional  shell  domain,  replace  nonlinear  differential  equations 
Eqs.  (6)  by  a set  of  two  nonlinear  algebraic  finite-difference  equations,  and 
solve  the  resulting  set  of  equations  by  the  nonlinear  relaxation  method. 

Details  of  the  nonlinear  relaxation  technique  will  not  be  discussed  here; 
description  of  the  technique  is  given  in  references  . Simply  put,  the 

nonlinear  relaxation  technique  offers  a method  of  systematically  reducing  the 
errors  at  each  nodal  point  for  each  algebraic  equation  to  some  acceptable  level. 

Rather  than  solving  6ill  two  differential  equations  simultaneously,  we  elect 
to  solve  first  the  w system  of  difference  equations  in  an  iterative  manner, 
and  subsequently  the  u system  of 'equations . We  revert  back  to  the  w system 
again  and  the  u and  so  on,  until  the  percentage  change  of  the  displacement 
at  all  nodal  points  is  always  less  than  0.0001  percent  on  the  (absolute) 
average . 

The  second  time  derivative  of  w in  Eq.  (6a)  is  approximated  by  the 
Houbolt's  third-order  backwards  difference  expression 

w(x,t)  = (l/S^)  |2w(x,t)  - 5w(x,t  - <5) 

+ W(x,T  - 26)  - w(x,T  - 36)j|  (9) 

where  6 = At  is  the  equal  time  increment. 

The  accuracy  of  this  representation  is  of  order  6^.  Special  attention  is 
devoted  to  the  first  few  time  steps  where  Eq.(9)  cannot  be  applied  directly. 
Before  giving  the  expression  f6r  these  first  few  time  steps,  we  note  that 
the  initial  conditions  are  of  the  form 


w(x,0)  = 0 


w(x,0)  = 0 


(10) 


i 


T - 


From  Eq.(lO),  expression  of  Eq.  ( 9)  for  the  first  few  time  steps  can 
readily  be  obtained  as  (Ref. 

1)  T = 0,  w{x,0)  = 0 

2)  T = 6;  since  w(x,0)  = 0,  we  have  w(x,-6)  = w(x,6) 
and  hence  w(x,6)  = (2/6^)  w(x,6) 

3)  T = 26;  w(x,26)  = (2/6^)  ^w(x,26)  - 3w(x,6)j 


k)  T > 26,  Eq.(9)  can  be  applied  directly. 


In  our  numerical  computation,  the  number  of  nodal  points  are  selected 
such  that  a subsequent  increase  in  nodal  points  does  not  significantly  affect 
the  magnitude  of  the  static  buckling  load.  With  this  consideration,  lU  nodal 
points  for  X = 5 and  up  to  22  points  for  X = 10  are  adopted;  A Poisson's 
ratio  of  1/3  is  also  used.  A time  step  must  be  selected  very  judiciously. 

A good  selection  is  such  that  the  results  are  within  a desired  accuracy,  but 
not  too  small  in  light  of  computer  time  considerations.  A reasonable  com- 
promise of  equal  time  increments  of  6 = 0.10  is  used;  this  same  time  increment 


is  also  selected  in  Ref 


■ ■ 
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DYNAMIC  BUCKLING  CRITERION 


Criteria  for  dynamic  axisymraetric  buckling  of  the  shallow  spherical 
shell  are  not  as  well  defined  as  for  static  buckling,  and  require  an  eval- 
uation of  the  transient  response  of  the  shell  for  various  load  levels. 


The  criterion  adopted  most  widely  (Refs.  10, 11,12, iHj  ) is  based  on 
plots  of  the  peak  nondimensional  average  displacement  in  time  history,  A max, 
of  the  shell  structure  versus  the  amplitude  of  the  load  where  A is  the 


average  displacement  and  has  been  defined 
follows : 


“d 


r a ra 

A = / rWdr  / / rzdr 


in  dimensionless  form  as 


(11) 


0 o 

The  numerator  is  the  volume  generated  by  the  shell  deformation  and  the 
denominator  is  the  constant  volume  under  the  cap. 

There  is  a load  range  where  a sharp  jump  in  peak  average  displacement 
occurs  for  a small  change  in  load  amplitude.  The  inflection  point  of  the 
load  deflection  curve  in  this  range  is  regarded  as  the  buckling  load;  this 
procedure  is  called  the  "inflection  point  method" . When  the  buckling  load 
is  defined  in  this  manner,  it  is  also  referred  to  as  the  dynamic  snap- 
through  load. 

The  method  is  demonstrated  by  a typical  example  shown  in  Fig.  2 . In 
this  figure,  A - t curves  are  displayed  for  several  uniform  load  parameters  p 
(this  example  is  related  to  a clamped  spherical  cap  of  X = 5 under  step 
uniform  pressure  of  infinite  duration).  A plot  of  vs.  p based  on 

information  given  in  Fig.  2 is  shown  in  Fig.  3 . It  is  clear  from  this 
figure  that  there  is  a sharp  jump  at  p = 0.U6,  and  according  to  the  criterion 
defined  above,  this  value  is  taken  as  buckling  pressure  for  this  particular 
spherical  cap  geometry  suid  load  situation. 
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CHECKS  WITH  EXISTING  SOLUTIONS  FOR  PERFECT  CAPS 

The  general  program  developed  for  buckling  of  spherical  caps  with  initial 
imperfections  [-]  is  modified  to  account  for  the  simpler  axisymraetric  spher- 
ical cap  under  dynamic  loading.  In  order  to  test  the  validity  of  this  modified 
program,  computer  runs  are  made  to  determine  dynamic  buckling  loads  of  spherical 
caps  for  which  niimerical  results  exist. 

Fig.  2 represents  A - x curves  for  different  load  levels  for  the  axisym- 
raetric spherical  cap  of  X = 5 under  uniform  step  pressure  of  infinite  duration; 
the  Ajjjgjj  - p curve  associated  with  Fig.  2 is  given  in  Fig.  3 . According  to 
the  dynamic  buckling  criterion  described  earlier,  dynamic  buckling  loads  for 
X = 5 under  step  loading  is  taken  as  p^^  = 0.U6  (where  is  denoted  as  the 
dynamic  buckling  or  snap-throxigh  load). 

Computer  runs  also  are  made  for  X = 7.5  and  10;  the  A - x curves  of  these 
two  geometric  parameters  are  shown  in  Figs,  h and  5 , respectively.  A 
summary  of  these  values  along  with  the  results  obtained  elsewhere  are  tabu- 
lated in  Table  1 for  comparison.  It  is  found  from  this  Table  that  the  present 
values  are  in  very  good  agreement  with  solutions  presented  in  Refs.  1^10-12, luj. 

Table  1 Dynamic  Buckling  Pressure  for  Various  Values  of  Spherical  Cap 
Parameters 


X 

5 

7.5 

10 

Present 

.1*6 

.1*1* 

.1*9 

Ref.  10 

.1*9 

.5 

.1*2 

Ref.  11 

.1*5 

.1*1* 

.37 

Ref.  12 

CO 

.51 

.50 

Ref.  ll* 

CD 

.5 

.1*3 
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The  present  solution  is  carried  out  using  time  steps  6 = 0.10  with 
lU  to  22  meridional  finite  difference  stations  for  X = 5 to  10.  The  length 
of  response  calculation  time  t carried  out  in  computer  nins  is  varied  from 
case  to  case  with  the  criterion  that,  in  the  neighborhood  of  the  buckling  load, 
T is  sufficiently  large  enough  to  allow  A - x ciirves  to  fully  develop.  As  a 
result  we  obtain  more  accurate  buckling  loauis,  and  save  computer  time  when 
load  levels  are  other  than  in  the  neighborhood  of  the  buckling  pressure. 

More  discussion  in  this  regard  will  be  provided  in  the  next  section. 

It  is  Interesting  to  note  from  Table  1 that  the  magnitude  of  dynamic 
buckling  loads  is  not  sensitive  to  shell  geometry  X.  This  finding  suggests 
that  the  dynamic  buckling  load  obtained  for  spherical  caps  represents  the 
value  for  the  complete  spherical  shell.  However,  this  observation  does  not 
apply  to  the  static  buckling  situation,  since  results  obtained  in  this  study 
show  that  static  buckling  loads  are  0.65,  1.02  and  0.85  for  X = 5,  7.5  and 
10,  respectively. 
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DYNAMIC  BUCKLING  OF  SPHERICAL  CAPS  WITH  INITIAL  IMPERFECTIONS 

With  the  validity  of  the  computer  program  established  for  dynamic  buckling 
analysis  as  demonstrated  in  the  previous  section,  we  now  proceed  to  inves- 
tigate the  effect  of  initial  imperfections  on  the  dynamic  buckling  capacity 
of  spherical  caps.  We  shall  also  examine  how  different  types  of  loadings  with 
various  time  duration  affect  the  dynamic  buckling  loads. 

Two  types  of  dynamic  loadings  are  considered  herein,  namely,  uniform  step 
pressure  of  infinite  duration  and  a right  triangular  pulse  with  various  time 
duration.  Solutions  obtained  in  this  study  will  give  some  indications  on  the 
most  severe  type  of  dynamic  loading  situations.  Specific  results  on  these 
two  types  of  lo8uiing  situations  are  discussed  separately  in  this  section. 

The  axisymmetric  initial  imperfection  adopted  in  this  study  is  of  the 
dimple  type  which  was  also  used  by  Koga  and  Hoff  ^2^.  This  type  of  imper- 
fection is  expressed  mathematically  as 

Wj  = (Wo/h)  (1  - x2)3  (12) 

where  W^  is  the  maximim  imperfection  which  occurs  at  the  shell  apex.  Selection 
of  this  expression  is,  in  fact,  quite  arbitrary;  however,  it  does  provide  an 
adequate  description  for  actual  shells  since  the  important  parameter  is  the 
maximum  eccentricity  and  not  the  imperfection  shape  function. 

Since  dynamic  buckling  pressure  is  not  sensitive  to  spherical  cap  geo- 
metric parameter  (X)  as  csin  be  seen  from  Table  1,  we  select  X * 5 for  all 
analyses  in  this  section. 


J 


r 
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Step  Loadins 


Let  us  first  consider  buckling  anstlysis  of  the  axisymmetric  spherical  cap 
(X  = 5)  with  initial  imperfections  xmder  a uniform  step  loading  of  infinite 
duration.  Some  typical  A - p curves  associated  with  various  load  levels  are 
given  in  Fig.  2 . From  this  figxire,  it  is  noted  that  a different  response 
calcualtion  time  t has  been  utilized  for  different  load  levels.  Although  the 
reason  for  using  this  strategy  was  mentioned  briefly  in  the  last  section,  a 
more  detailed  discussion  in  this  regard  is  given  here.  As  is  apparent  frcan 
Fig.  2 the  peak  value  of  A occurs  after  several  cycles  of  oscillation  for 
the  load  level  p above  the  critical  value  of  p for  snap-through  (we  shall  call 
this  value  p^p  the  dynamic  buckling  or  dynamic  snap-through  load).  It  is 
also  observed  that  for  the  load  levels  above  Per*  larger  the  value  of 
p,  the  smaller  the  response  calculation  time  taken  for  A to  reach  its  peak 
value.  Clearly,  a flexible  response  calculation  time  is  desired  to  cope  with 
veurious  load  situations. 

The  strategy  initiated  here  stems  from  the  fact  that  in  dynamic  analysis 
with  significant  nonlinearity,  time  consiming  iteration  schemes  and  a huge 
number  of  time  steps  are  usually  Involved.  This  strategy  perhaps  can  best 
be  accomplished  through  the  execution  of  the  computer  program  in  an  interactive 


mode. 


By  inserting  some  simple  statements  in  the  program,  execution  of 


the  program  will  pause  after  a fixed  period  of  response  calculation  time  and 
the  user  can  decide  to  continue  or  terminate  the  execution  depending  on 
whether  A has  reached  its  peak  vsilue.  On  the  other  hand,  for  a load  range  very 
close  to  Pqp  , sufficient  response  calcualtion  time  should  be  given  to  Insure 
that  the  shell  structure  does  not  eventually  buckle. 

This  approach  should  predict  more  accurate  buckling  loads  and  will  also 
save  considerable  computing  time.  Huang  carried  his  solution  for  X * 5 

out  to  T = lU  using  a very  fine  time  increment  6 = 0.02  and  (UX  + 1)  stations. 
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i 

llj  obcainecl  their  solucior  up  to  t = 60  for  1 = 5 
and  T = 120  for  higher  values  of  X by  using  6 = 0.10  and  20  stations.  Ball  ^ 

lU^  carried  their  results  out  to  t = 50  for  X = 5,  and  t = 120  f 

for  X = 8 and  11  stations  with  6 = 0.05.  As  mentioned  earlier,  in  the  ’ 

present  calculations,  we  employ  different  response  calculation  times  for  * 

different  load  situations;  6 = 0.10  and  ih  to  22  stations  for  X = 5 to  10 
are  selected.  The  convergence  criterion  adopted  is  that  the  average  abso- 
lute percentage  change  of  displacement  functions  is  less  than  0.0001.  i 

With  adoption  of  this  strategy  for  the  selection  of  response  calculation  | 

time,  a series  of  runs  are  made  for  an  axisymmetric  spherical  cap  (X  = 5)  ' ^ 

with  different  magnitudes  of  initial  imperfections.  Fig.  6 represents  the  _ 

^raax  “ ^ curves  for  this  shell  geometry  with  imperfections  W^/h  ranging  from  i 

I 1 

0.1  to  1.0.  Aslo  Included  in  this  figure  for  comparison  purposes  is  the 
result  for  the  perfect  shell  which  has  already  been  displayed  in  Fig.  3 . 

According  to  the  dynamic  buckling  criterion  described  earlier  in  this  paper, 
the  dynamic  buckling  loads  (p^p)  are  found  to  be  0.U6,  0.39,  0.28  and  0.185 
for  imperfections  of  WQ/h  = 0,  0.1,  0.5  and  1.0,  respectively. 

Sudden  jump  in  A near  the  critical  load  area  is  very  obvious  for 
W^/h  = 0,  0.1  and  0.5.  For  the  case  of  W^/h  = 1.0,  the  -p  curve  does 

not  exhibit  such  a drastic  change  in  displacement  and  the  inflection  point 
technique  is  used  to  determine  the  magnitude  of  the  buckling  load.  Suffice 
it  to  say  that  the  sudden  Jump  phenomenon  may  be  viewed  as  a particular  case 
of  the  situation  whose  dynamic  buckling  load  must  be  determined  by  the  inflec- 
tion point  method.  Therefore,  the  inflection  point  method  can  be  regarded  ■ 

as  a more  general  means  to  obtain  dynamic  buckling  loads. 

A plot  of  P(,p  vs.  W^/h  based  on  the  results  displayed  in  Fig.  6 is 
shown  in  Fig.  T . From  this  figure  it  is  obvious  that  initial  imperfections 
do  indeed  have  the  effect  of  reducing  dynamic  buckling  capacity  of  spherical  caps. 
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when  comparing  with  dynamic  buckling  loads,  static  response  yields  higher 
buckling  values  for  W^/h  in  between  0 and  0.7**,  but,  on  the  other  hwd,  has 
a lower  buckling  veilue  when  W^/h  is  greater  than  O.jk,  Evldentally,  when 
imperfections  exist,  the  static  response  reduces  shell  buckling  capacity  at 
a faster  rate  than  that  by  dynamic  response.  In  other  words,  static  response 
of  axisymmetric  spherical  caps  is  more  sensitive  to  initial  Imperfections 
than  dynamic  response.  This  finding  is  significant  in  shell  stinicture  design 
if  initial  imperfections  are  tedten  into  consideration. 
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Triangular  Pulse 

Having  analyzed  the  dynamic  buckling  response  of  the  spherical  cap  with 
an  axisymmetric  imperfection  Eind  a spatially  and  temporally  constant  pressure 
field,  we  turn  next  to  the  same  problem  except  for  the  pressure  pulse  which 
is  now  taken  as  triangular  rather  than  constant  with  time.  The  right  tri- 
angualr  pulse  shown  in  the  subset  of  Fig.  8 is  taken  as  representative  of 
this  type  of  loading.  At  dimensionless  time  t*  the  pressure  has  a zero  value. 

Except  for  some  minor  modifications  made  to  deal  with  the  time-varying 
nature  of  applied  loading,  the  same  computer  program  and  strategy  described 
earlier  are  utilized  for  this  analysis. 

Four  triangle  pulses  with  duration  x**  = 0.5^^,  2.176,  10.88  and  21.76 
are  considered.  The  same  spherical  cap  geometry  (X  = 5)  und  magnitude  of 
imperfections  are  retained.  Fig.  8 presents  the  results  of  the  axisymmetrical 
spherical  cap  under  triangle  pulse  of  x*  = 0.5^^  and  with  imperfections 
W^/h  = 0,  0.1,  0.5  and  1.  According  to  the  dynamic  buckling  criterion  adopted 
in  this  paper,  buckling  loads  are  discerned  as  5.05,  ^.8,  3.95,  and  3.1  for 

= 0,  0.1,  0.5,  and  1,  respectively.  Figs.  9 » 10  , and  11  display 
the  results  associated  with  triangular  pulse  durations  x*  = 2.176,  10.88  and 
21.76,  respectively. 

A summary  of  numerical  results  shown  in  these  figures  is  given  in  Table 
2 for  comparison.  Also  included  in  this  table  are  the  values  for  static 
response  and  step  loading  of  infinite  duration.  In  addition  to  buckling 
loads.  Table  2 also  gives  for  most  situations  the  dimensionless  buckling  time 
and  finite  difference  node  where  plasticity  begins.  One  interesting  and 
important  observation  made  from  this  table  is  that,  for  the  cases  of  very 
short  triangular  pulses  (x^*  = 0.5^^  and  2.176),  the  spherical  cap  buckles 


i 


after  removal  of  the  pulses.  Evidently,  the  kinetic  energy  imparted  to  the 
cap  by  the  pressure  pulse  was  enough  to  propel  the  cap  into  a buckling  mode 


I] 
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1 

well  after  load  removal.  An  examination  of  L nodes  provided  in  this  table  j 

^ 1 

reveals  that  no  regular  pattern  is  developed  for  the  location  as  to  where  | 

i 

plasticity  may  begin  (Mises  criterion  utilized).  ! 

’< 

Fig.  13  exhibits  the  relationship  between  dynamic  buckling  load  p^^  > 

» * 

and  pulse  duration  t*  for  different  magnitudes  of  imperfections.  This  f 

figure  demonstrates  that  the  length  of  has  a very  significant  effect  on 
the  magnitude  of  the  dynamic  buckling  load.  For  the  same  W^/h,  this  figure 
shows  that  dynamic  buckling  load  for  = 0.544  is  about  9 to  11  times 
greater  than  that  for  = 21.76. 

Fig.  14  gives  the  curves  for  dynamic  buckling  load  versus  initial 
imperfection  for  various  triangular  pulse  durations.  Also  shown  for  com- 
parison purposes  is  the  result  of  step  loading  which  was  displayed  in  Fig.  7. 

As  can  be  seen  from  this  figure  that,  for  a given  magnitude  of  imperfection, 
the  shorter  the  triangular  pulse,  the  higher  the  buckling  load.  This  figure 
also  reveals  that  the  limiting  case  of  triangular  pulse  of  various  time 
durations  is  the  step  loading  of  infinite  duration.  It  is  important  to 
point  out  that  this  limiting  case  represents  the  most  severe  loading  situa- 
tion for  d3niamic  buckling  analysis,  since  it  yields  the  lowest  values  of 
dynamic  buckling  loads. 
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DISCUSSION  AND  CONCLUSIONS 

An  investigation  has  been  carried  out  for  dynamic  buckling  analysis  of 
ajci symmetrical  spherical  caps  with  initial  Imperfections.  Two  types  of  load- 
ing are  considered,  namely,  step  pressure  of  infinite  duration  and  triangular 
pulses  of  various  time  durations.  Governing  equations  given  in  Refs, 
for  buckling  euisLlysis  of  spherlceLL  caps  with  asymmetric  imperfections  are 
modified  to  account  for  the  axisymmetric  nature  of  the  problem  and  dynamic 
effects.  Central  finite  differences  and  Houbolt's  backwards  difference  scheme 
Eire  used  to  replace  spatial  and  temporal  derivatives,  respectively.  The 
nonlinear  relaxation  method  ^l6,17^  is  then  employed  to  solve  the  resulting 
algebraic  equations. 

Unlike  static  buckling  aneilyses,  dynamic  buckling  loads  are  deter- 


H 


H 


mined  from  a displaced  volume-pressure  - p)  curve  by  the  inflection  point 


H- 


method  I 111.  Dynamic  snap-through  occurs  when  there  is  a sharp  Jump  in  the 
Voc  - P H-  This  indicates  that  the  maximum  amplitude  of  oscil- 

lation of  the  shell  has  a dlscontlnnous  change. 

To  verify  possible  validity  of  the  present  computer  program,  a series 
of  runs  are  made  to  determine  dynamic  buckling  loads  for  axisymmetrlcal 
sphericeJ.  caps  of  X = 5»  7.5,  and  10  under  uniform  step  loading  of  infinite 
time  duration.  Current  results  along  with  previous  findings  ^10-ll»^  are 
tabulated  in  Table  1 for  comparison.  Good  agreement  among  these  solutions 
is  observed;  it  may  be  concluded  that  the  inflection  point  method  has  been 
generally  accepted  as  a dynamic  buckling  criterion. 

Since  dynamic  buckling  loads  are  not  sensitive  to  shell  geometry  as  can 
be  seen  from  Table  1,  we  select  X = 5 for  all  cuialyses  in  this  paper.  A 
dimple  type  of  imperfection  is  also  adopted  in  this  study  which  provides  a 
very  adequate  description  of  the  local  nature  of  initial  imperfections  in 
spherical  shells. 


» 
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Results  of  the  axisymmetric  spherical  cap  (X  = 5)  under  step  loading  are 
obtained  for  majcimum  imperfections  W^/h  = 0.1,  0.5  and  1.0.  These  results 
are  displayed  in  Fig.  6 along  with  that  of  the  perfect  shell.  Buckling  loads 
obtained  from  Fig.  6 by  the  inflection  point  method  are  presented  in  Fig.  7 . 
Also  displayed  in  Fig.  7 • for  comparison  purposes  are  the  solutions  from 
static  buckling  analysis.  It  is  observed  from  this  figure  that,  for  W^/h  less 
than  about  0.7^,  the  static  buckling  load  is  higher  than  the  dynamic  buckling 
load.  However,  when  W^/h  is  greater  than  0.7^,  the  above  situation  is 
reversed.  This  result  demonstrates  that  the  static  buckling  load  decreases 
more  rapidly  than  the  dynamic  buckling  load  because  of  imperfections.  Moreover, 
the  static  response  is  more  sensitive  to  initial  imperfections  than  the  dynamic 
response. 

It  is  important  to  point  out  that  imperfections  do  indeed  have  the  effect 
of  reducing  buckling  capacity  for  both  static  and  dynamic  responses,  although 
they  are  affected  in  different  manners. 

Calculations  are  also  carried  0’'t  for  the  same  shell  geometry  with  the 
same  magnitude  of  imperfections  under  right  triangular  pulses.  Four  different 
pulse  durations  are  considered,  namely,  t*  = 0.5^^,  2.1T6,  10.88  and  21.76. 
Results  of  A - p for  these  four  cases  are  shown  in  Figs.  8-11  , respectively. 
Dynamic  buckling  loads  obtained  from  these  figures  are  displayed  against  pulse 
duration  t*  in  Fig.  13  for  various  magnitudes  of  imperfection  and  against 
initial  imperfection  in  Fig.  lU  for  different  pulse  durations.  Also  displayed 
in  Fig.  lU  is  the  solution  for  the  step  loading.  From  Fig.  13  , it  is  found 
that  pulse  duration  has  a very  significant  impact  on  the  magnitude  of  dynamic 
buckling  load.  For  example,  for  the  same  imperfection,  dynamic  buckling  load 
for  T*  = 0.5^^  is  about  10  times  higher  than  that  for  x*  = 21.76.  This  obser- 
vation can  eilso  be  seen  from  Fig.  lU  . In  fact.  Fig.  l4  shows  that  the  limiting 
case  of  the  triangular  pulse  is  the  step  loading  of  infinite  duration,  and  that 
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the  step  loading  represents  the  most  severe  loading  situation  for  dynamic 
buckling  analysis,  since  it  yields  the  lowest  values  of  buckling  loads. 

A summary  of  results  displayed  in  Fig.  8-11  along  with  those  of  static 
and  dynamic  step  loadings  of  infinite  duration  is  given  in  Table  2 for  com- 
parison purposes.  In  addition  to  buckling  loads.  Table  2 also  provides 
information  for  the  buckling  time  and  finite  difference  node  where  plas- 
ticity begins.  Two  observations  are  madet  first,  for  the  cases  of  very 
short  triangular  pulses  ( t*  = 0.544  and  2.176),  the  spherical  cap  buckles 
after  removal  of  the  pulses;  secondly,  no  regular  pattern  is  developed  for 
the  location  as  to  where  plasticity  may  begin. 

As  a final  remark,  the  flexible  response  calculation  time  option  set 
up  in  the  computer  program  enables  users  to  continue  or  terminate  the  exe- 
cution. This  option  not  only  provides  a possibility  of  obtaining  more  ac- 
curate dynamic  buckling  loads,  but  also  saves  considerable  computing  time 
for  the  cases  where  additional  execution  appears  to  be  unnecessary. 
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Fig.  2 Dynamic  response  of  a clamped  spherical  cap  ( X = 5)  under  uniform  step  loading  of  infinite  duration 
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Fig.  3 Dynamic  buckling  load  for  a clamped  spherical  cap  ( X = 5)  under  uniform  step  loading  of 
infinite  duration. 


Fig.  4 Dynamic  response  for  a clamped  spherical  cap  ( X « 7.5)  under  uniform  step  loading  of  infinite  duration. 


Fig.  6 Dynamic  buckling  loads  for  a clamped  spherical  cap  ( X = 5)  with  axisymroetric  initial  imperfections 
and  under  uniform  step  pressure  of  infinite  duration. 
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and  under  triangular  pulse  ( 't  = .544). 


2.4 


a 


Fig.  13  Dynamic  buckling  load  vs  triangular  pulse  duration  for  an  axi symmetric 
spherical  cap  (^=  5)  with  initial  imperfection. 
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